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I. INTRODUCTION 



Despite their obvious physical difference, from a mathematical (or purely formal) point of 
view, magnetism and superconductivity are closely linked phenomena. Coherent charge and 
pairing fluctuations can be described in an abstract "charge space" in terms of precession of 
pseudospins operators, first introduced by Anderson, and behaving as spin-1/2 operators^!. 
The charge part of the fluctuations couples to the square of an external electric field through 
Raman matrix elements, opening the possibility to manipulate pseudospins by light pulses in 
a way that is analogous to nuclear magnetic resonance (NMR) and electron-paramagnetic- 
resonance (EPR) experiments for real spins 3 . 

In our experiment, a polarized ultrafast laser pulse excites the superconductor^' through 
the Impulsive Stimulated Raman Scattering (ISRS) effect 8 . The coherent oscillations of the 
Cooper pairs condensate are detected via delayed supercontinuum pulses and enable a new 
technique, Coherent Charge Fluctuation Spectroscopy (CCFS), to distinguish the electronic 
excitations that couple to the superconducting quasiparticles. This is of pivotal importance 
for cuprates, as the applicability of conventional pairing theories-, based on retarded in- 
teractions between electrons mediated by low energy glue bosons, has been doubtecP^ 
and a completely different framework has been proposed involving non-retarded interactions 
associated with high-energy electronic scales^. 

These high-temporal resolution (< 50 fs) experiments were performed in two optimally 
doped (T c = 40 K) La 2 _ a; Sr. E Cu04 (LSCO) single crystals (x = 0.15) with different orienta- 
tions, see RefP^and Appendix for details. A polarized 1.55 eV laser pulse with a duration of 
45 fs and an absorbed fluence around 300 /x J / cm 2 induces both dipole (linear in the electric 
field) and Raman (quadratic in the electric field) allowed excitations, the latter being the 
main focus of this work. 

II. EXPERIMENTS 

We chose different experimental geometries for exploiting the Raman selection rules {Ap- 
pendix) for excitation and detection and obtain information on different final states. In the 
first geometry, the pump electric field is parallel to the Cu-0 bond giving access to Raman 
excitations with Ai g + B\ g symmetry, while the probe pulse electric field is directed towards 
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the c-axis, allowing us to detect only A\ g symmetry excitations. Then, using the same 
pumping geometry, we probed the excited system along [100] and [010], which respectively 
give access to A\ g + Bi g and A\ g — Bi g excitations. Performing the difference between the 
two orientations allows us to extract only B\ g excitations. Finally, we used the pumping and 
probing fields on the diagonal direction giving access to A\ g + B2 9 Raman excitations. The 
dynamics of all these excitations is then probed by broad-band ultrafast reflectivity, whose 
overall time-energy dependence is displayed in Fig. [TJA-C. 

III. RESULTS AND DISCUSSION 

The transient reflectivity is dominated by a large abrupt amplitude change followed by 
a relaxation; this is a consequence of high energy particle-hole (p-h) excitations produced 
by the dipole allowed absorption of the pump photons. Furthermore, in both orientations 
the transient reflectivity changes sign throughout the spectra in correspondence to specific 
electronic transitions. These changes reflect the transfer of spectral weight among the dif- 
ferent absorption bands produced by the high energy p-h excitations 6 . The number of p-h 
excitations involved is estimated in the Appendix to be less than 10~ 2 per Cu atom (Fig. [4]). 

All geometries present coherent oscillations of Raman excitations on top of the dipole p-h 
excitations relaxation. In Ai g symmetry, an ultrafast oscillation with a period of 145 fs and 
a long coherence time (1.45 ps) is visible at all wavelengths (Fig. [l]^4)- 

The Fourier transform analysis of A± g symmetry data is presented as an inset in Fig. [Tp. 
A sharp peak at 28 meV is visible, corresponding to the out-of-plane La Ai g mode of LSCO^. 
Such coherent fully symmetric modes have already been observed in high-T c superconduc- 
tor*!™. 

Instead, in both A lg + B 2g and B\ g symmetries (Fig. [TjE and difference in Fig. [l]F), 
slower and damped (around 300 fs coherence time) oscillations are clearly observed below 
T c . The Fourier analysis of these time-resolved profiles is presented in Fig. [2]A-i?. In Ai g +B2 9 
symmetry, for a probing wavelength of 2.45 eV, the temperature dependence of the Fourier- 
transform signal shows an obvious peak at 18 meV which vanishes above T c (Fig. [2]S). In B lg 
symmetry, the broad peak appears at energies around 24 meV when the sample temperature 
is lower than T c . Increasing the pump fluence to 2 mJ/cm 2 , no such peak could be observed 
below T c (Fig.||A). 
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FIG. 1: Transient broadband reflectivity data at 10 K in A\ g [pump || [100], probe || [001], (A, 
D)], A lg + B 2g [pump || [110], probe || [110], {B, E)} and at 24 K A lg + B lg [pump || [100], probe || 
[100], (C, F)] geometries (specified in tetragonal axis). The extracted profiles are shown in panels 
(C, F) for selected probe energies. Panel (E) presents the reflectivity oscillations by substracting 
the background on the profile, and in panel (F) we show the difference between A\ g + B\ g and 
A\ g — Bi g profiles, which is proportional to the B\ g signal. The absorbed pump fluence is around 
300 u,J/cm 2 . 

In Fig. [2]C-.D we display the THz spectra obtained in the superconducting state, and 
compare them with the spontaneous Raman response (data taken frorrP^), which is well un- 
derstood in terms of the excitation of two Bogoliubov quasiparticles^. The good agreement 
between them allows us to identify the strongly temperature dependent part of the oscilla- 
tions as Raman charge fluctuations of the superconducting condensate. This agreement is 
expected from simple theoretical considerations (Appendix) for electronic ISRS which show 
that any excitation that is Raman active in a collinear configuration of incoming and out- 
going photon electric field is also accessible in a pump-probe experiment. Our experiment 
detect remnant of superconductivity at fluences of the same order but larger than previ- 
ously reportecP^. We attribute that difference to the much higher sensitivity of our optical 
measurement to superconductivity and its bulk character. Probably superconductivity is 
indeed quenched on the first layers of the sample and becomes invisible to surface probes 
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like photoemissiori^. 
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FIG. 2: Fourier transform spectra obtained at different temperatures and excitation nuences in B\ g 
[A, 1.91 eV probing energy) and A\ g + B2 g (B, 2.45 eV probing energy) geometries; Panels (C) and 
(D) show the comparison between transient reflectivity data and Raman measurements, in the su- 
perconducting phase. The spontaneous Raman spectra are the difference between superconducting 
and normal phases, showing only the charge fluctuation peaks. The insets show schematically: in 
panel (C), the angle (j> along the Fermi surface and shows the regions in momentum space excited 
in Big symmetry; panel (D), idem for £>2g symmetry. 

The temporal evolution of the coherent phonon oscillation in A\ g geometry is presented 
in Fig. [l]D, with its extrapolation down to zero-time delay (as defined irP^), allowing us 
to establish its cosine waveform which is typical of a displacive (resonant) mechanism of 
excitation! 19 * 20 *. The electronic transitions induced by 1.55 eV photons occurs between the 
ground state of the material and higher energy electronic states. At this energy, a peak in 
the optical absorption is observed in LSCO, distinct from the Cu-0 charge transfer^, which 
had been attributed to charge ordering in the form of stripes^. Thus, the cosine wave form 
indicates that, not unexpectedly, the charge modulations are strongly coupled to the A\ g La 
phonon. 

In the A\ g geometry, the presence of the strong coherent phonon disturbs the real-time 
observation of the superconducting condensate. Instead, in both A lg + B 2g and B lg symme- 
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tries, the fluctuations of the superconducting quasiparticles are clearly observable and start 
at zero time delay, allowing the determination of a sine waveform (Fig. [IjE'-F). This indicates 
that contrary to the A\ g phonon case, the triggering mechanism is Impulsive Stimulated Ra- 
man Scattering^^, meaning that the p-h excitations absorbed by the pump pulse are not 
directly coupled to the superconducting quasiparticles. We show below that an analysis of 
the probe energy dependence leads to the same conclusion. 

For the Raman allowed excitations, the effect of the pump on the electrons can be de- 
scribed by a time-dependent impulsive potential quadratic in the electric fielcP^ (Appendix). 
Its effect on the superconducting quasiparticles is easily understandable using Anderson's 
pseudospin formalism 1 . We define^ the pseudospin operators cr k behaving as spin-1/2 op- 
erators and such that an up (down) pseudospin means that the states k f and — k 4- are 
simultaneously empty (occupied), while a sidewise pseudospin means particle hole mixing 
which is the hallmark of quasiparticles participating in the pairing. The reduced BCS 
Hamiltonian corresponding to a <i-wave superconductor in the presence of a time- dependent 
potential can be written in term of the cr k operators as 

H BCS = -J2K + sw(t)]. ( r ]i (i) 

k 

where we introduced the fictitious "magnetic fields" b k = (A k ,0,£ k ) an d <5b k (t) = 
(0, 0, v£ (t)), with A k the momentum dependent superconducting order parameter at equi- 
librium, £ k the quasiparticle band energy (in the absence of superconductivity) measured 
from the Fermi level and fjf (£) the Raman part of the pump induced time dependent ex- 
ternal potential (Appendix). In the ground state, for v£(t) = 0, each pseudospin aligns 
parallel to b°. In the <i-wave state the pseudospins form a texture (Fig. [3]/l) with sideway 
pseudospins close to the Fermi level (£ = 0) and away from the nodal point indicating regions 
participating in the pairing. 

The pseudospins obey the equations of motion 

h d ^ = -2K + 6b k (t)] x<r k . (2) 

implying that for t > the pseudospins precess around the equilibrium direction with an 
angular velocity 2|b k |//i, |b k | = y^JT A£ being the BCS quasiparticle energy (see Fig. ^3 
and Fig. [5|. Since <5b k is in the z direction, only pseudospins having a significant component 
at equilibrium in the x-y plane respond to the Raman impulsive field, automatically selecting 
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FIG. 3: Pseudospin textures coding the BCS wave function in momentum space. (A) Ground 
state texture: pseudospins are labeled by the distance in energy from the Fermi level £ = £k and 
the angle <f> along the Fermi surface (inset of Fig. 2 C); (B) Small red arrows: amplitude of the 
impulsive field 5h^ applied at t = in B<i g symmetry. Long arrows: texture snapshots (amplitudes 
exaggerated for clarity) immediately after the excitation (yellow), at 25 fs (green) and at 50 fs 
(blue); (C) Theoretical charge fluctuation: solid dots corresponds to the snapshots of panel (B). 
The open dots are the experimental change in reflectivity after the high energy ph background has 
been subtracted; (D) Probe energy dependence of the Fourier transformed A\ g + Big fluctuation. 

the quasiparticles participating in the pairing, thus in Bi g symmetry only pseudospins which 
are close to the Fermi level and neither in the nodes nor in the antinodes have a significant 
time dependence. The z projection of the pseudospin, in Fig. [3]B, encodes the contribution 
to the total charge fluctuation shown in panel C. At t = 50 fs (blue) the pseudospins at 
£ = and close to <fi = 7r/4±7r/8 are close to their maximum negative amplitude in the z 
direction, corresponding to the first peak in the B 2g charge fluctuation. 

We also show in Fig. [3]C a comparison between the experimental and theoretical con- 
densate oscillation in Big. Interestingly the experiment show a quite long coherence time 
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compared to theory. Details of the computations are given in the Appendix. 

The transient optical properties of the system in the presence of a fluctuation of symmetry 
X = Ai g , Bi g , B 2g are governed by the changes in the dielectric function tensor 



and dx/d(Nx) is the conventional Raman tensor (Appendix) . We see that the same Raman 
tensor appears in the generation of the pulse by ISRS and in the subsequent probe. In 



the fluctuating quasiparticles will contribute to dxl (Nx) allowing to detect excitations 
participating in the pairing. We point out that CCFS is not restricted to reflectivity, and 
other techniques like spontaneous Raman scattering can be used as a probe allowing to 
test also excitations of different symmetry. In this case a different matrix element will be 
involved in the probe. 

The oscillation of the superconducting condensate is most clearly visible in the A\ g + B2 g 
configuration. For this reason, we perform the spectral analysis in this configuration. The 
probe-energy dependence of the A% g + B 2g fluctuation in the frequency domain is presented 
in Fig. |3p. The superconducting fluctuations clearly resonate at an energy of 2.6 eV, cor- 
responding to the Cu-0 charge transfer energy of the parent compound which coincides 
with the Hubbard U energy of a one-band description^. Remarkably, even though there is 
substantial absorption below the charge transfer band in our samples, the superconducting 
quasiparticles appear to be decoupled from the excitations in the measured range. This is 
fully consistent with our finding above that the Ai g + B 2g fluctuations have a sine waveform. 

The correct framework to understand superconductivity in cuprates has been subject 
of an intense debate^ 1 11 1 25 1 . One possibility is that the role of phonons in the traditional 
mechanism is replaced by a different low energy bosonic excitation like damped magnons 
which act as a glue allowing electrons to pairing^E3 i n this scenario, superconductivity 
can be understood by the traditional frameworlP where retardation plays an important role. 
Anderson^ has argued that there is no such a low energy glue and that proximity to the 
Mott phase is an essential ingredient. The relevant time scale of the interactions producing 
the pairing is the inverse of the Hubbard energy U ^ 2eV. Therefore, the interaction can be 
considered instantaneous for practical purposes. Our results are consistent with a coupling of 
the superconducting quasiparticles with excitations at 2.6 eV which are a clear fingerprint of 




(3) 



analogy with lattice ISRSP^, only excitations having an interaction matrix elements with 
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"Mottness" in the superconducting state. A negligible coupling in the rest of the measured 
energy window (1.6 eV < hu < 3.2eV) is observed, but we cannot exclude that other 
electronic excitations outside our probing range are also coupled to superconductivity and 
even dominant. Numerical computations support a coupling to the Mott scale^, although 
with a strong contribution from the low energy region. 

The experimental breakthrough leading to the conclusion that phonons act as the pairing 
glue in conventional superconductors was the isotope effectpSES. The key feature of this 
experiment is its high specificity, since only the frequency of one potential glue excitation 
is affected and its impact on superconductivity evaluated. CCFS has a high degree of 
specificity in a reverse form, only paired electrons are affected and its impact on different 
excitations assessed. The presented results form a benchmark for time-resolved experiments 
in cuprates and shed new light on the nature of the pairing interactions. 

In a more general perspective the NMR/EPR analogy encoded in Eqs. ([[]), ([2| allows to 
borrow concepts^ like the relaxation times T\ and T 2 * . T 2 is defined by the decay of the 
charge fluctuations, which is dominated by the inhomogeneity of the pseudomagnetic field 
in momentum space. Therefore, our experiment opens appealing perspectives to typical 
NMR/EPR-like techniques such as coherent control of the superconducting wave-function 
by sequence of pulses. These tools can be generally applied to different materials including 
heavy fermions and iron-based superconductors. 
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Appendix A: Experiment 

Two single crystals of optimally doped LSCO were cut, polished and oriented via X- 
ray diffraction to obtain two surfaces containing respectively the a and b or the a and c 
crystallographic directions. 
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Pump-probe reflectivity was performed with a 1 kHz Ti:Saph amplified fs laser whose 
output was used to generate a white light continuum probe (1.5-3.2 eV) and a monochromatic 
(1.55 eV) pump beam, capable of fluences up to 10 mJ/cm 5 ^. The pulse duration is 45 fs 
which puts a lower limitFto the frequency of the excitations that can be excited with ISRS, 
u> > 1/(45/*). 

The reflectivity spectrum and a reference signal were synchronously detected by two 
identical spectrometers. 

Appendix B: Effet of particle-hole excitations on the optical properties 

Here, we estimate the changes produced in the optical properties of the sample due to the 
creation of dipole allowed particle-hole (p-h) excitations by the pump pulse. Within a mean- 
field picture, the optical conductivity or the charge susceptibility x °f the sample changes 
due to: i) change of the initial state and final state occupation and ii) modification of the 
electronic structure due to the out of equilibrium distribution. We present a computation 
of the second effect. We assume that the main outcome of the p-h excitations is to change 
the balance between the Cu (rid) and O (n p ) hole occupation numbers putting the charge 
transfer (CT) hole density Tier = 2n p — out of equilibrium. To take into account Cu 
d x 2_ y 2 and O p XjV orbitals, the system is described by a three-band Hubbard model: 

H = S ^e l n i + Yl tijiClCjc + h.c.) 

i <ij>cr 

+ 2J Uiftifhn, + 2J Uijfiihj. 

i <ij> 

Here Cj a , Ci a are creation and destruction operators for holes on lattice site i with spin a 
and we defined the occupation number operators, n ia = Cj a Ci a , hi = h^ + h^. The on-site 
energies are £j = e p (e<f) and the on-site repulsions, Ui = U p (Ud) for i in an O (Cu) site; and 
the hopping matrix elements, = ±t p d(±t pp ) for Cu-0 (O-O) nearest neighbor sites, with 
the sign depending on the relative phase of the Wannier orbitals involved. For the inter-site 
interaction, we keep only the nearest neighbor repulsion Uij = U p d between Cu and O. We 
define A = e p — We take the same parameters as in RefP, except for A = 3eV. This 
ensures that the first maximum in the charge transfer band in the optical conductivity of 
the parent compound is at 2.16 eV, in good agreement with the experiments^, instead of a 
higher value obtained with the first principle parameters of RefP. 
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To account for the strong correlations on Cu, we treat the d orbital using the Gutzwiller 
approximation and the other interactions using the Hartree-Fock approximation. This leads 
to a single particle Hamiltonian with renormalized energies, which at equilibrium depend self- 
consistently on the charge distribution. Our strategy is to compute the optical conductivity 
at mean-field level in the ground state and in the presence of the out of equilibrium charge 
distribution. We also neglect any asymmetry in this out-of-equilibrium distribution. In 
other words we assume that the out-of-equilibrium distribution has A\ g symmetry. 

The leading effect of the out-of-equilibrium charge distribution is to renormalize the 
effective energies at mean-field level. Since this happen through total charges that are 
integrals of the distribution, this effect will be rather insensitive to its precise form. The 
charge is put out of equilibrium by adding a constraint implemented through a Lagrange 
multiplier. Thereafter, the optical conductivity and the charge susceptibility x are obtained 
at mean-field level, and the derivative respect to the out of equilibrium charge is evaluated 
numerically. For the doped case, the ground-state is assumed to consist of stripes, as in 
RefP, which gives a good overall agreement with experiments for the optical conductivity. 

Figure [4] shows the result of the computation for the transient imaginary part of the 
in-plane dielectric function obtained as 



with 5ncT(t) = ncrit) — n° CT and n° CT being the equilibrium charge transfer density. 

Given the simplifications in the computation, we obtain a fair overall agreement with the 
experiment adjusting only the parameter #n<7r(500/,s) = 8 • to fit the intensity. This 
gives the estimate of the number of transferred holes per Cu. 

The number of excited electrons per unit cell may be estimated independently in terms 
of experimental parameters using the following formula: 



where A is the absorption coefficient and l s the penetration depth, both at the pumping 
wavelength, V the unit cell volume, F the pumping fluence in J / m 2 and AE the transition 
energy. For a fluence of 300 fiJ/cm 2 , we obtain a value of n = 4 • 10 -3 per Cu atom, in fair 
agreement with the theoretical calculation of excited electrons Sricr- 

Sncrit) is positive, which is consistent with the fact that in the ground state there are far 
more holes on Cu than on the O, so the laser pulse will tend to decrease and increase n p . 



de xx (u),t) = -Ait 



dn CT 



(u)Sncr(t)- 
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FIG. 4: Comparison between the experimental (red thick line) and theoretical (blue thin line) 
transient imaginary part of the in-plane component of the dielectric function tensor £2 at 500 fs 
time delay for doping x = 0.15. For the theory we used <5ncr(500/s) = 8 • 10~ 3 and three band 
parameters as specified in the text. 

The transient dielectric function has a strong energy dependence due to transfer of spectral 
weight among the different absorption bands which leads to regions of positive and negative 
transient reflectivity as shown in Fig. 1 A-B. 

The profile of e 2 can be understood as due to a softening of the charge transfer edge due to 
a downward renormalization of the mean-field charge transfer energy A. The leading effect 
is due to the nearest neighbor repulsion 5 A ~ —2U p dSncT- The structures appear at higher 
energy than in the experiment probably because our mean-field approach underestimates 
the softening of the CT edge with doping. 

The same strategy can be used to compute the electronic Raman tensors defined in the 
main text and in the previous section. 

Appendix C: Coherent generation of Cooper-pair condensate charge fluctuations 

In this section, we discuss the computation of the Raman charge fluctuation in the su- 
perconducting state. Within a one-band description of electrons close to the Fermi surface 
we consider uniform (i.e. zero momentum) charge fluctuations described by the operator^, 
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kcr 

Here f£ l9 = [cos(k x a)+cos(k y a)]/2 1 /jf l9 = [cos(k x a)—cos(k y a)]/2 : f^ 2g = sin(k x a) sm(k y a) , 
n ka = c kcr c kcr is the occupation operator for the state with wave vector k and spin a and c' ka 
(ck a ) are creation (destruction) operators for electrons. 



Generalizing the arguments of Ref.^21 for lattice ISRS to electronic ISRS, we write the 
Hamiltonian of the system in the presence of the pump pulse as, H = Hbcs + Hr where 
Hbcs describes the low energy quasiparticles in a BCS ground state and 

x 

is the perturbation due to the pump laser. 

In the semiclassical approximation and for ul much larger than the frequency of the 
fluctuations, 

vx{t) = ~l E{t) ■ w E(t) - (S1) 

Here x is the charge susceptibility, dxif^h) /d{Nx) is the conventional second rank Raman 
tensor for electronic scattering with symmetry X and incident frequency cuzp^, E is the time 
dependent electric field of the pump wave and X runs over the allowed symmetries. 



In spontaneous Raman scattering the operator has the same form as in Eq. (SI) but 
with electric fields which can be change independently on the left and on the right of the 
Raman tensor. Thus the selection rules for ISRS and spontaneous Raman are quite similar. 
Disregarding orthorombicity and using the symmetry properties of the Raman tensor in 
LSCO 18 (Z?4/i group) for a pump pulse polarized in the plane the Raman operator reads 



H R 



° K -' -\{e x f + {e y f]N Alg 



[d(N Alg ) 

(S2) 



with e a versor in the direction of the electric field (E = Ee) and axes parallel to the 
CuO bond direction. We see that for a pulse in the [100] direction, A\ g + Bi g symmetries 
are excited [Fig. 1 A and C\ while if the while if the electric field is aligned along [110], 
A\ g + B 2g symmetries are excited [Fig. 1 B] . Selection rules for the detection can be derived 
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FIG. 5: Snapshots of the textures coding the BCS wave function in momentum space. Small red 
arrows: amplitude of the impulsive field <5bk applied at t = in B\ g symmetry. Long arrows: 
texture snapshots (amplitudes exaggerated for clarity) immediately after the excitation (yellow), 
at 25 fs (green) and at 50 fs (blue). 

analogously from Eq. (3). In the case of Fig. 1A the probe field is in the [001] direction and 
only the A\ g component of the A\ g + B\ g fluctuation can be seen. 
For symmetry X the fluctuations at zero temperature are given by, 

-t 



(Nx)(t) 



I dt'J^^Mt - t')]\(0\N x \u)\ 2 v x (t') (S3) 

J — oo 



with the sum running over a complete set of states. If the pump pulse of width r and height 



e is approximated by a Dirac delta function {vxif) = 5(t)re), Eq. (S3) yields a sine like wave 
form with shape determined by the Fourier transform of the conventional Raman scattering 
line shape for symmetry X^. 

We compute the dynamics of the charge fluctuation in a d-wave BCS state using pseu- 
dospins operators defined as 

°k = (cktC-14 + h.c), 
14 



i&i = (cktC-14 - h.c), 
a\. = 1 - n kt - 71-14. 

Using the definitions in the main text we have that each pseudospin is subject to a time 
dependent potential fjf (£) = vx(f)f£- 

We first linearize the equation of motion [Eq. (2)] for the pesudospins^^ in terms of the 
time dependent fluctuations <5cr k (t) = cr k (t) — cr k with cr k the equilibrium texture shown in 
Fig. 3 A The linearized equation reads, (2) 

= -2b k x <5<r k - 2<5b k x < 

We define the versor e_i_ = y x cr k in the direction perpendicular to b k and the y direction 
and pseudospin projections <5cr k = <5er k e_i_ + <5<r k y. Axes are defined on Fig. 3B and on Fig. [5} 
For simplicity, we set the temperature to zero and neglect collective effects and Coulomb 
interactions which can be easily incorporated in the random phase approximation^^. 

Solving the pseudospin equation of motion for a time dependent impulsive potential 
applied at t — 0, v*(t) = 8{t)ref* one obtains, 

5^ = -2sin(2^)^re/^ 

K = 2cos(2 J E; k t)^re/ k Y 

which describe the pseudospins precessing around b k . 

Figures [5] and 3 B show the evolution of the texture in Bi g and B 2g symmetry respec- 
tively. Notice that only quasiparticles participating in the pairing have a significant time 
dependence and these are further restricted by the symmetry of the impulsive potential. 

The coherent charge fluctuation is determined by the component of the oscillations in the 
z direction: 2 

N x (t) = 2j2^m(^^) re t > 0. (S4) 

k V k / 

The sum runs over all the Brillouin zone. The squared factor above selects only the paired 
quasiparticles further restricted by the symmetry function. This shows again the selectivity 
of CCFS. 



15 



Eq. (S4) has been evaluated numerically using a <i-wave gap function Ak = A [cos(fc x .a) — 
cos(k y a)}/2 with A = 20meV and the one band parameterization of the electronic structure 
of LSCO given in Ref.^. The result for B 2g symmetry is shown in Fig. 3 C. 
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